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Abstract. The family Hylodidae is composed of 46 species distributed in three genera: Crossodactylus (13 species), Hylodes 
(26 species), and Megaelosia (seven species). Although the monophyly of the Hylodidae is supported by previous molecu- 
lar phylogenetic inferences, the monophyly of Megaelosia and/or Hylodes has been questioned. Crossodactylus and Hylodes 
share the plesiomorphic diploid chromosomal number 2n = 26, whereas in the species of Megaelosia as karyotyped up to 
now, diploid numbers have ranged from 28 to 32. Here, we expand to six (of seven) the number of species of Megaelosia 
sampled in a phylogenetic approach based on mitochondrial and nuclear DNA sequences and describe the 26-chromo- 
some karyotype of M. goeldii for the first time. Our results provide an improved perspective on the relationships among 
the frogs described previously in the genus Megaelosia, and we propose a new taxonomic arrangement, in which Mega- 
elosia is a monotypic genus, with a new genus being described to accommodate the remaining species allocated previously 
to Megaelosia. We also confirm that 2n = 26 is the plesiomorphic diploid number in the Hylodidae, and conclude that an 
increased diploid chromosome number, together with two morphological traits and one ethological characteristic, are po- 


tential synapomorphies of the new genus. 
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Introduction 


The family Hylodidae is currently composed of 46 spe- 
cies distributed in three genera: Crossodactylus DUMERIL 
& BIBRON, 1841 (13 species), Hylodes FITZINGER, 1826 (26 
species), and Megaelosia MIRANDA-RIBEIRO, 1923, with 
seven species (FROST 2021, SEGALLA et al. 2021). The ge- 
nus Megaelosia is currently composed of seven species 
- Megaelosia apuana POMBAL, PRADO & CANEDO, 2003, 
M. bocainensis GIARETTA, BOKERMANN & HADDAD, 1993, 
M. boticariana GIARETTA & AGUIAR-JR., 1998, M. goeldii 
(BAUMANN, 1912), M. jordanensis (HEYER, 1983), M. lutzae 


IZECKSOHN & GOUVEA, 1987, and M. massarti (DE WITTE, 
1930) (FROST 2021). In addition to the overall microen- 
demism of most Megaelosia species, they typically occur 
at low population densities (TOLEDO et al. 2014) and are 
known only from their type localities or extremely limited 
areas of occurrence (GIARETTA et al. 1993, SANTOS et al. 
2011, Muscar et al. 2020a), located primarily within the 
Serra da Mantiqueira and Serra do Mar mountain rang- 
es, in southeastern Brazil [in the states of São Paulo (SP), 
Rio de Janeiro (RJ), Espirito Santo (ES), and Minas Gerais 
(MG)]. By contrast, a number of populations of M. goel- 
dii have been found at multiple localities of the Brazilian 
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states of Sao Paulo and Rio de Janeiro (HADDAD et al. 2013, 
FROST 2021). 

The phylogenetic inferences published on the Hylodidae 
up to now have included a reduced number of species, and 
have not supported the monophyly of either Hylodes (Py- 
RON & WIENS 2011) or Megaelosia (GRANT et al. 2017) due 
to the relationships of M. goeldii. In the phylogenetic tree 
provided by PyRON & WIENS (2011), the genus Hylodes was 
recovered as paraphyletic with respect to M. goeldii, which 
was the only Megaelosia species in the analysis. In the phy- 
logenetic proposal by Grant et al. (2017), which included 
three Megaelosia species, M. goeldii was closely related to 
the species of the genus Hylodes, rendering Megaelosia a 
paraphyletic genus. More recently, the phylogenetic analy- 
sis by Dusots et al. (2021) also included M. goeldii as the 
only representative species of its genus, and proposed that 
Megaelosia is a junior synonym of Hylodes. This synonymi- 
zation was reversed by SEGALLA et al. (2021), who revali- 
dated the genus Megaelosia until further and more con- 
clusive analyses were available. Given this scenario, a phy- 
logenetic analysis based on a more representative sample 
of the genus Megaelosia will be needed to better evaluate 
the interspecific and intergeneric relationships within the 
Hylodidae. 

Hylodes and Crossodactylus are cytogenetically similar. 
Except for Hylodes nasus, which has a diploid number of 
2n = 24 (BOGART 1991), the karyotypes of the Hylodes and 
Crossodactylus species all have the same diploid number 
(2n = 26) and a highly similar chromosome morphology 
(BEGAK 1968, BRUM-ZORRILA & SAEZ 1968, BOGART 1970, 
DENARO 1972, DE Lucca & JIM 1974, AGUIAR-JR. et al. 
2004, AMARO 2005). In the genus Megaelosia, on the other 
hand, diploid numbers range from 2n = 28 to 2n = 32 chro- 
mosomes (MELO et al. 1995, GIARETTA & AGUIAR-JR. 1998, 
Rosa et al. 2003). That is, while M. lutzae has a diploid 
number of 2n = 32 (ROSA ET AL. 2003), and M. boticaria- 
na is 2n = 30 (GIARETTA & AGUIAR-JR. 1998, Rosa et al. 
2003), M. massarti has 2n = 28 or 29, due to the presence 
of B chromosomes, that is, Bı or B2 (MELO et al. 1995, 
Rosa et al. 2003). The karyotype of M. goeldii, which is 
the species that has represented its genus in many of the 
published phylogenies, was unknown prior to the present 
study. 

In order to advance our knowledge on the genus Mega- 
elosia, in the present study, we (i) describe the karyotype of 
M. goeldii and (ii) verify the monophyly of this genus based 
on mitochondrial and nuclear DNA sequences. Finally, as 
a result of our analyses we describe a new genus, to avoid 
paraphyly in the genus Megaelosia. 


Material and methods 
Specimens analyzed 


For our cytogenetic analysis, we sampled four Megaelosia 
goeldii tadpoles, which had been collected at the type local- 
ity (Table 1). The samples used for the phylogenetic infer- 
ences were obtained from tadpoles or adults of M. goeldii, 


M. apuana, M. boticariana, M. massarti, M. jordanensis, 
and M. lutzae, which were collected, as far as possible, at 
their respective type localities (see Table 1, Fig. 1). We also 
included individuals collected from Pindamonhangaba 
(Fazenda Votorantin Celulose) in São Paulo state (SMRP 
81.3, SMRP 81.6, and SMRP 81.8 - see Supplementary Sup- 
plementary Table S1), which are the same samples used by 
Rosa et al. (2003), who identified the species as M. lutzae. 
As Pindamonhangaba is about 100 km (in a straight- 
line) from Itatiaia National Park, RJ, the type locality of 
M. lutzae (IZECKSOHN & GOUVEA 1987), and said authors 
provided no information on the criteria they used to iden- 
tify the species as M. lutzae, we conservatively treated these 
specimens as Megaelosia sp. 

The ingroup in the phylogenetic analyses also includ- 
ed samples of two adults of Hylodes asper from Santo An- 
dré, two adults of H. japi from Jundiai, and one adult of 
H. phyllodes from Santo André, all in the state of São Paulo 
(SP); one adult of H. heyeri from Morretes, in the state of 
Parana (PR); and two adults of Crossodactylus gaudichau- 
dii from Rio de Janeiro (RJ). For information on museum 
catalogue numbers, collecting localities, and the GenBank 
accession numbers of the nucleotide sequences, see Sup- 
plementary Table S1. 

All our samples of cells and tissues were deposited in 
the Shirlei Maria Recco-Pimentel (SMRP) collection of 
the Laboratory of Chromosomal Studies (LabEsC) and the 
Luis Felipe Toledo tissue collection (TLFT), both in the 
Institute of Biology at Campinas State University (UNI- 
CAMP), Brazil. The collection of specimens was author- 
ized by the Brazilian Institute for the Environment and Re- 
newable Natural Resources (IBAMA #02027.011919-95), 
through the Chico Mendes Institute for the Conservation 
of Biodiversity (ICMBio/SISBIO #17914-2/21619-1, #52352- 
2 and #75770-1), and were registered in the National Sys- 
tem for the Management of Genetic Patrimony and As- 
sociated Traditional Knowledge (SisGen #ACD362E). For 
tissue extraction, the specimens were euthanized through 
the application of anaesthetic (10% Lidocaine) to the skin, 
following the recommendations of the Herpetological Ani- 
mal Care and Use Committee (HACC) of the American 
Society of Ichthyologists and Herpetologists (available at 
http://www.asih.org). 

To complement our matrix, we downloaded GenBank 
sequences of M. goeldii, M. jordanensis, M. boticariana, 
H. phyllodes, H. nasus, H. meridionalis, H. caete, H. sazi- 
mai, H. perere, H. ornatus, H. amnicola, H. lateristrigatus, 
Crossodactylus schmidti, and C. caramaschii. The outgroup 
included species representative of the family Alsodidae, 
which has been inferred to be the sister taxon to Hylodidae 
(PYRON & WIENS 2011, GRANT et al. 2017; for details, see 
Supplementary Table S1). 


Description of the karyotype 


Mitotic metaphases were obtained from cell suspensions 
of the intestinal epithelium of four specimens of M. goeldii 
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Table 1. Species of Megaelosia analysed in the present study and their respective sampling localities in Brazil. The type locality of each 


species, where pertinent, is indicated by an asterisk. 


Species Specimens analysed Locality 
Megaelosia goeldii 4 tadpoles Serra dos Órgãos National Park, Teresópolis, RJ* 
Megaelosia apuana 3 tadpoles Caparaó National Park, Alto Caparaó, MG 


Megaelosia boticariana 3 tadpoles 


Gruta Funda Park, Atibaia, SP* and Sao Francisco Xavier, São José dos Campos, SP 


Megaelosia massarti 2 adults Paranapiacaba, Santo André, SP* 
Megaelosia jordanensis 1 tadpole Campos do Jordão, SP* 
Megaelosia lutzae 5 tadpoles Itatiaia National Park, Itatiaia, RJ* 
Megaelosia sp. 3 tadpoles 


Fazenda Votorantin Celulose, Pindamonhangaba, SP 


treated previously with colchicine (KING & ROFE 1976, 
with the modifications of GATTO et al. 2018). The chro- 
mosomes were stained with Giemsa (10%) and C-banded 
(KiNG 1980), and then silver-stained using the Ag-NOR 
method (HowELL & BLACK 1980) to detect the Nucleo- 
lus Organizer Regions (NORs). We used 2n as the somatic 
chromosome number, and FN as the fundamental number 
or number of chromosome arms, following WHITE (1954). 
The metaphasic chromosomes were photographed under 
an Olympus BX-60 microscope (Tokyo, Japan), edited in 
Adobe Photoshop, and classified according to GREEN & 
SESSIONS (1991). 


4PW 


21°S 


23°S 


25°S 


Extraction of the DNA 


Genomic DNA was extracted from liver or muscle tis- 
sue, previously stored in ethanol at -80°C. The tissue was 
lysed in TNES (50 mM Tris pH 7.5, 400 mM NaCl, 20 mM 
EDTA, and 0.5% SDS) supplemented with proteinase K 
(100 pg/mL) at 56°C for approximately 4 hours. After lysis, 
the samples were treated with NaCl to a final concentration 
of ~1.7 M. The DNA was precipitated in isopropyl alcohol, 
washed in ethanol (70%), and rehydrated in TE (10 mM 
Tris-HCl, 1 mM EDTA, pH 8). For quality control, and 
to quantify the genomic DNA, the samples were electro- 
phoresed in 0.8% agarose gel and analyzed by NanoDrop 
(Thermo Fisher Scientific). 


Hypsometry 
Om 


2890 m 


Figure 1. Schematic map showing the geographic distribution of the Megaelosia samples analysed in the present study. Red star = 
Megaelosia goeldii (Teresópolis, RJ), with its geographic distribution shaded in red; blue triangle = Megaelosia apuana (Alto do Caparaó, 
MG); blue circle (filled) = Megaelosia jordanensis (Campos do Jordão, SP); blue circle (open) = Megaelosia sp. (Pindamonhangaba, 
SP); blue rhombus = Megaelosia boticariana (Atibaia, São José dos Campos and Caçapava, SP); blue pentagon = Megaelosia massarti 
(Piracicaba, SP), and blue square = Megaelosia lutzae. Filled symbols indicate type localities. 
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Sequencing of the mitochondrial and nuclear DNA 


To generate data for the phylogenetic analysis and ge- 
netic distance estimates, a mitochondrial fragment that 
included the 12S and 16S rRNA genes and the tRNA-Val 
gene was amplified by PCR using the primer pairs MVZ 
59 (GRAYBEAL 1997) and Titus I (TrTus 1992) (or MVZs50; 
GRAYBEAL 1997), and 12SL13 (FELLER & HEDGES 1998) (or 
HıKF — CARVALHO ET AL. 2020) and 16Sbr (PALUMBI et 
al. 2002). In some cases, we used the primers 12SKf and 
HiKR (CARVALHO et al. 2020) to amplify a fragment of 
the 12S rRNA gene. A segment of the nuclear RAG-1 gene 
was also amplified by PCR, using the primers RAG-1F and 
RAG-1R (FAIVoVICH et al. 2005). The amplified products 
were electrophoresed in 1% agarose gel and then purified 
using the GFX PCR and Gel Band DNA Purification kit 
(GE Healthcare) according to manufacturer’s instructions. 
The samples were sequenced using the BigDye Terminator 
kit (Applied Biosystems). In addition to the primers cited 
above, we used the internal primers 16SL2a, 16H10 (HEDG- 
ES 1994), and 16Sar (PALUMBI et al. 2002) for sequencing. 

The products of the sequencing reaction were puri- 
fied by precipitation in 80% ethanol and centrifuged at 
1,200 rpm for 30 minutes, and were then washed in 70% 
ethanol and centrifuged for 10 minutes. Once dried, the 
products were resuspended in loading dye (Blue-Dextran- 
EDTA/Formamide, 1:5), denatured for 3 minutes at 94°C 
and then transferred to an automatic sequencer. The se- 
quences were edited using the Bioedit software, available at 
http://www.jwbrown.mbio.ncsu.edu/BioEdit/bioedit.html 
(HALL 1999). 


Analysis of the DNA sequences 


The mitochondrial matrix (2,337 bp) was concatenated 
with that of the Ragi fragment (428 bp). The sequences 
were aligned using MAFFT v7 (KaToH et al. 2019; https:// 
mafft.cbrc.jp/alignment/server/), following the G-INS-i 
strategy, and the resulting matrix had 2,765 bp and 68 ter- 
minals. 

The phylogenetic analyses were based on Maximum 
Parsimony (MP) and Bayesian Inference (BI) approach- 
es. The MP analyses were conducted in TNT v1.5 (GoLo- 
BOFF et al. 2016) using the new technology search option 
(including sectorial searches, ratchet, tree drifting, and 
tree fusing). Gaps were considered as a fifth state, and 
node support was evaluated by bootstrap analysis with 
1,000 replicates. The BI was run in BEAST v2.6 (BouUCK- 
AERT et al. 2019), based on the GTR model of DNA evo- 
lution, inferred by MrModeltest v2.3 (NYLANDER 2004). 
A relaxed clock log normal approach and the Yule proc- 
ess were applied. The final Markov Chain was run twice 
for 100 million generations. Each tree was sampled every 
10,000 generations and the first 25% of the trees were dis- 
carded as burn-ins. Tracer v1.5 was used to check the con- 
vergence of the chains (RAMBAUT & DRUMMOND 2009). 
The trees were summarized with TreeAnnotator v6.0 and 


visualized in FigTree v1.4.2 (RAMBAUT & DRUMMOND 
2019). 

Confidence intervals of the support values were cal- 
culated following HUELSENBECK & HILLis (1993) and 
FELSENSTEIN (2004), and the nodes with a posterior prob- 
ability (PP) of > 0.95 and bootstrap (BS) values of at least 
75% were considered to be sufficiently resolved, whereas 
those with a PP of 0.90-0.95 and BS values of 50-75% were 
considered to indicate a tendency, and those with a PP of 
less than 0.90 or BS below 50% were considered to be un- 
resolved. 

The genetic distances (p-distance) between species were 
calculated using a 1,402-bp fragment of the mitochondrial 
16S gene (delimited by the primers 12SL13 and 16SBr) and 
a 538-bp fragment of the 16S gene (the 16Sar-16SBr frag- 
ment). All ambiguous positions were removed from each 
pair (pairwise deletion option) and the analyses were run 
in MEGA X (Kumar et al. 2018). 


Inference of ancestral diploid numbers 


To infer the ancestral chromosome number in Megaelosia 
we used the ChromEvol server (GLICK & MAYROSE 2014, 
RICE & MAYROSE 2020) and the Bayesian Inference topol- 
ogy from this work as the input. To find the model that 
best fits our data, we first performed a test using all models 
available on ChromEvol and selected the model with the 
lowest Akaike Information Criterion (AIC) value, which 
was “BASE NUM”. A second test was performed with op- 
timized parameters in the best model. The haploid chro- 
mosome numbers of M. goeldii (from this work), M. boti- 
cariana, M. lutzae, M. massarti, three Crossodactylus spe- 
cies, and five Hylodes species (as previously cited) were 
plotted into the analysis. Unknown chromosome numbers 


«> 


were coded “x”. 
Nomenclatural acts 


The electronic edition of this article conforms to the re- 
quirements of the amended International Code of Zoo- 
logical Nomenclature, and hence the new names contained 
herein are available under that Code from the electronic 
edition of this article. This published work and the nomen- 
clatural acts it contains have been registered in ZooBank, 
the online registration system of the ICZN. The LSID (Life 
Science Identifier) for this publication is: urn:lsid:zoobank. 
org:pub: 7FF6FDDC-4869-4C87-9FC3-26359042BEE6. 


Results 
Megaelosia goeldii karyotype description 


The karyotype of M. goeldii has a diploid number of 2n = 
26 and a fundamental number of FN = 52, being com- 
posed of one metacentric chromosome pair (pair 4) and 
12 submetacentric pairs (pairs 1-3 and 5-13) (Fig. 2a). In 
the Giemsa-stained metaphases, secondary constrictions 
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Table 2. Uncorrected p-distances (in percentages) between M. goeldii and the Phantasmarana species based on the 1,402-bp (below 
diagonal) and 538-bp fragments (above diagonal) of the mitochondrial 16S gene. The grey cells in the diagonal show intraspecific 


distances (larger fragment/smaller fragment). 


Taxon 1 2 3 4 5 6 7 
1. Megaelosia goeldii 0.1/0.2 8.2 73 7.2 7.5 6.9 6.7 
2. Phantasmarana apuana 13.4 0/0 1.7 1.6 1.6 3.2 3.5 
3. Phantasmarana sp. 13.1 1.9 0.2/0.2 0.2 1.2 2.6 2.7 
4. Phantasmarana jordanensis 13.0 1.7 0.3 0.3/0.2 1.1 2.5 2.6 
5. Phantasmarana lutzae 13.1 2.1 1.7 1.6 0.3/0.4 2.8 2.8 
6. Phantasmarana massarti 12.6 4.5 4.0 3.9 4.0 0.2/0 0.9 
7. Phantasmarana boticariana 12.5 4.3 3.9 3.8 4.0 2.3 0.9/0.1 


were observed in the pericentromeric region of the long 
arm of the homologs of pair 3 (Fig. 2a), which coincided 
with the NORs detected by silver-staining (Fig. 2b). Size 
heteromorphy of the homologous NORs was observed in 
two of the four individuals analyzed. The C-bands were de- 
tected in the centromeric region of all the chromosomes 
and pericentromerically in the long arm of the homologs 
of pair 3 (coinciding with the NORs) and in the short arm 
of pair 12 (Fig. 2c). 


Phylogenetic inferences on the Hylodidae 


In both the BI and the MP analyses, the internal relation- 
ships in the monophyletic family Hylodidae recovered 
paraphyly of Megaelosia in relation to Hylodes (Fig. 3), 
with M. goeldii being the sister taxon of the clade com- 
posed of the 12 Hylodes species included in our dataset. 
The M. goeldii+Hylodes clade had moderate support (pos- 
terior probability = 0.9, bootstrap = 68%). All the remain- 
ing species of Megaelosia analyzed here (i.e., M. massarti, 
M. boticariana, M. apuana, M. jordanensis, M. lutzae, and 
Megaelosia sp.) were recovered in a highly supported clade 
(posterior probability = 1.0, bootstrap = 100%), which was 
the sister group of the M. goeldii+Hylodes clade (Fig. 3). 
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This topology indicates the need for a taxonomic rear- 
rangement to avoid paraphyly in Megaelosia and to this 
end, we propose a new genus, which is described below. 

The new genus presented well-supported internal rela- 
tionships in our phylogenetic analyses, with M. massarti as 
the sister species of M. boticariana, and M. jordanensis as 
the sister species of M. apuana and M. lutzae (Fig. 3). 

The specimens SMRP 81.3, SMRP 81.6, and SMRP 81.8, 
identified a priori as Megaelosia sp., and analyzed previ- 
ously by Rosa et al. (2003), were nested within our sample 
of M. jordanensis, rather than the topotypes of M. lutzae, 
in both topologies (Fig. 3). The genetic distances between 
these three Megaelosia sp. individuals and the M. jordanen- 
sis specimens, including a M. jordanensis topotype (CFBH 
28578), estimated from a 1,402-bp fragment of the 16S 
rRNA gene, were very low (approximately 0.3%; Table 2). 
We thus conclude that the specimens collected by Rosa et 
al. (2003) are, in fact, individuals of M. jordanensis. 


Ancestral diploid numbers 
The ancestral chromosomal number analysis recovered 


the node that includes Hylodes+M. goeldii with 2n = 26 
(posterior probability = 0.98). The node that includes the 


n 


e 


" 


12 


Figure 2. Karyotype of Megaelosia goeldii. (a) Giemsa staining, (b) Chromosome pair 3 showing the nucleolus organizer regions 
(NORs) detected by silver-staining, (c) C-banding. Scale bar = 5 um. 
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other Megaelosia species (i.e., M. massarti, M. boticariana, 
M. apuana, M. jordanensis, M. lutzae, and Megaelosia sp.) 
was recovered with 2n = 28 as the most likely ancestral dip- 
loid number (posterior probability = 0.75) for these spe- 
cies. 


Eupsophus 2n= 28-30 


100 
= Alsodes 2n= 22-34 
100 
0.6: 
0.6: 
0.98, 
0.99/86 
0.99, 
100 
0.99; 
0.90 
0.58/99 


Phantasmarana gen. nov. 


ZooBank: LSID urn:lsid:zoobank.org:pub:E26B7159-0F73-4D09- 
AB4E-E39E24FoB23F 


Type species: Elosia massarti DE WITTE, 1930 


Crossodactylus gaudichaudii ZUEC 1 757 | 
Crossodactylus gaudichaudii ZUEC 175 
Crossodactylus trachystomus KY0071 — 
Crossodactylus schmidti AY843579 

Crossodactylus caramaschii KI961569 —— 
Crossodactylus werneri KU215901 
Crossodactylus werneri KU215900 
Hylodes nasus KJ961577 AE 
Hylodes asper ZUEC 11555 o w 
Hylodes asper ZUEC 11441 

Hylodes heyeri ZUEC 24650 

Hylodes meridionalis MF624225 

Hylodes meridionalis KT221614 

Hylodes lateristrigatus KM390) 

oo Hylodes caete KY627909 

Hylodes caete KY627908 


Hylodes phyllodes ZUEC 115: 
Hylodes phyllodes DQ50217i 
Hylodes phyllodes DQ283096 


Hylodes ornatus KJ961578 "Em 
Hylodes japi ZUEC 24643 

Hylodes japi ZUEC 24642 

Hylodes sazimai KJ961585 

Hylodes sazimai KJ961584 

Hylodes amnicola KJ961576 

Hylodes amnicola KJ961575 

Hylodes perere KJ961581 
Hylodes perere KJ961580 
Megaelosia goeldii SMRP 527. 
Megaelosia goeldii SMRP 527. 
Megaelosia goeldii SMRP 527.. 
Megaelosia goeldii DQ28307. 
Phantasmarana massarti ZUEC ‘ao, an 
Phantasmarana massarti ZUEC 1 PE O) 5 
Phantasmarana boticariana SMRP 10! 
Phantasmarana boticariana SMRP 109 
Phantasmarana boticariana KJ961586 
Phantasmarana boticariana ZUEC 24. 
Phantasmarana jordanensis CFBH 2857 
Phantasmarana jordanensis SMRP 81 
Phantasmarana jordanensis MF6242. 
Phantasmarana jordanensis SMRP 81 
Phantasmarana jordanensis SMRP 81.3 
Phantasmarana apuana SMRP 526.4 
Phantasmarana apuana TLFT 1981 
Phantasmarana apuana TLFT 1911 
Phantasmarana lutzae 555.2 
Phantasmarana lutzae 555.1 
Phantasmarana lutzae 555.4 
Phantasmarana lutzae 555.5 
Phantasmarana lutzae 555.3 


100 


Figure 3. Phylogenetic relationships of the family Hylodidae inferred by Bayesian analysis of the mitochondrial sequence H1 and 
the RAG-1 fragment. The numbers at the nodes are posterior probabilities (left) and bootstrap values in percentage (right), inferred 
from Bayesian and Maximum Parsimony analyses, respectively (values of less than 50% have been omitted). The asterisk indicates 
Phantasmarana sp., which we assigned to P. jordanensis. Diploid numbers (2n) were obtained in previous studies (cited in the text; 
BLOTTO et al. 2012, QUERCIA et al. 2020), except for M. goeldii, which is described in the present study. 
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Diagnosis. This genus was recognized as monophyletic 
in our phylogenetic analyses based on mitochondrial and 
nuclear DNA sequences and can also be diagnosed by the 
combination of the following phenotypic traits: (i) pres- 
ence of paired lateral white vocal sacs (vs. absent in Mega- 
elosia, completely translucent in Hylodes, and partially 
translucent in Crossodactylus: Fig. 4), (ii) presence of a vo- 
cal slit (vs. absent in Megaelosia), (iii) fold of fifth toe does 
not reach outer metatarsal tubercle (vs. reaching this tu- 
bercle in Megaelosia, not extending much beyond the base 
of the proximal subarticular tubercle in Crossodactylus and 
Hylodes), (iv) snout rounded or slightly protruding in lat- 
eral view (vs. clearly protruding in lateral view in Megaelo- 
sia), (v) finger discs more than 30.5 times the SVL (vs. less 
than 30.5 times in Megaelosia), (vi) scutes on finger discs 
absent or weakly developed (vs. well developed in Meg- 
aelosia), (vii) adult snout-vent length exceeding 6 cm (vs. 
less than 5 cm in Hylodes and Crossodactylus), (viii) ab- 
sence of advertisement call (vs. present in Crossodactylus 
and Hylodes, except for H. vanzolinii; but absent in Meg- 
aelosia), and (ix) diploid numbers with 28 or more chro- 
mosomes (vs. 26 in Megaelosia, Hylodes and Crossodacty- 
lus). Traits i, iii, and ix represent putative synapomorphies 
for this new genus. 


Content. Six species. Phantasmarana apuana (POMBAL, 
PRADO & CANEDO, 2003) comb. nov.; P bocainensis (G1A- 
RETTA, BOKERMANN & HADDAD, 1993) comb. nov.; P. boti- 
cariana (GIARETTA & AGUIAR-JR., 1998) comb. nov.; P. jor- 
danensis (HEYER, 1983) comb. nov.; P. lutzae (IZECKSOHN 
& GOUVEA, 1987) comb. nov.; and P massarti (DE WITTE, 
1930) comb. nov. 


Sister taxon. The clade including Hylodes and Megaelosia. 


Description. Adult frogs of this genus are large, with SVL 
ranging from 6 to 12 cm in adult males and females. Tad- 
poles are giant, reaching 13 cm or more in total length (Sup- 
plementary Table S2 and Appendix S1) and may remain in 
the larval stage for more than one year (AUGUSTO-ALVES et 
al. unpublished data). Adults and larvae are likely noctur- 
nal, inhabit forest streams, and adults are generalist, feed- 
ing on invertebrates, fish, and other amphibians (PELOSO & 
PAVAN 2007, ALVES & TOLEDO 2017, AUGUSTO-ALVES et al. 
2018). Although they lack an advertisement call, a mysteri- 
ous type of call has been reported in P. boticariana (Mus- 
CAT et al. 20208), and adult males may use their contrast- 
ing white vocal sac for visual signalling in intraspecific 
nocturnal communication (AUGUSTO-ALVES et al. 2018). 


Figure 4. Adult male Crossodactylus caramaschii (frame extracted from a video by D. PERRELLA: ZUEC-VID 800) highlighting its 
partially translucent vocal sac (A); Adult male Hylodes phyllodes (photograph by L. F. ToLepo) and Hylodes asper (photograph by M. 
Martins) with completely translucent vocal sacs (B-C); preserved male P. boticariana (D: photograph by L. F. ToLepo: ZUEC 11843) 
and live male P. apuana (E: frame extracted from video by G. Aucusto-ALvEs: ZUEC-VID 402; AucusTo-ALvEs et al. 2018) with 
contrasting white vocal sacs. 
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Etymology. The name Phantasmarana is derived from the 
combination of Latin words phantasma (= phantom) + 
rana (frog), meaning “phantom frogs” in reference to the 
extraordinary rarity of frogs of this genus in the wild (only 
few herpetologists have seen these frogs in their natural 
habitats), the lack of an advertisement call, and the fact that 
some enigmatic sounds have been reported. The gender is 
feminine. 


Distribution. The genus Phantasmarana occurs in forest- 
ed areas of the Atlantic rainforest biome, ranging from the 
state of São Paulo (southern limit) to the states of Minas 
Gerais and Espirito Santo (northern limit). Endemic to 
Brazil. 


Megaelosia MIRANDA-RIBEIRO, 1923 
Type species: Hylodes goeldii BAUMANN, 1912 


Diagnosis. This genus is monotypic and can be diagnosed 
by the combination of the following traits: (i) males lacking 
vocal sacs and vocal slits (vs. present in Hylodes and Phan- 
tasmarana); (ii) fold of fifth toe reaching outer metatarsal 
tubercle (vs. not reaching this tubercle in Phantasmarana); 
(iii) snout subacuminate in dorsal view and protruding in 
profile (vs. rounded or slightly protruding in lateral view in 
Phantasmarana)); (iv) tips of fingers and toes with a pair of 
scutes free from subunguis (vs. absent or weakly developed 
in Phantasmarana); (v) finger discs less than 30.5 times 
the SVL (vs. more than 30.5 times in Phantasmarana); 
(vi) prootic not fused with the frontoparietals (vs. fused in 
Hylodes, unknown in Phantasmarana); (vii) alary process 
of the hyoid absent (vs. rudimentary in Hylodes, unknown 
in Phantasmarana); (viii) mesosternum composed of car- 
tilaginous plaques (vs. bony, broad, and bifurcated posteri- 
orly in Hylodes, unknown in Phantasmarana); (ix) absence 
of apical supplementary elements of the m. intermandibu- 
laris (vs. present in Hylodes and Phantasmarana)); (x) zygo- 
matic ramus larger than optic ramus (vs. absent in Hylodes, 
unknown in Phantasmarana); (xi) absence of advertise- 
ment call (vs. present in Hylodes spp., except in H. vanzoli- 
nii), and (xii) diploid number = 26 chromosomes (vs. 28 or 
more in Phantasmarana). Traits i and ix (lack of vocal slit, 
sac, and supplementary elements of the m. intermandibu- 
laris) are autapomorphies of this genus. 


Content. One species: Megaelosia goeldii (BAUMANN, 1912). 
Sister taxon. The genus Hylodes. 


Description. Adult frogs of this genus are large, with SVL 
varying between 8.2 and 9.5 cm SVL in both adult males 
and females (Supplementary Table S2). Tadpoles are large, 
one was measured, and found to be 11.7 cm in total length 
(NuIN 2003). Adults and larvae are thought to be diurnal 
(GIARETTA et al. 1993, NUIN 2003) and adults are general- 
ists, feeding on invertebrates and other amphibians (G1a- 


RETTA et al. 1993). Although there is no advertisement 
call, an unidentified type of call was reported in M. goeldii 
(Muscat et al. 2020b). 


Distribution. The single species of the genus occurs in for- 
ested areas of the Atlantic rainforest biome, from the north 
of the state of Sao Paulo (Serra da Bocaina) to central Rio 
de Janeiro state (Desengano State Park). Endemic to Brazil. 


Discussion 
Intrageneric relationships in the family Hylodidae 


The family Hylodidae now consists of four genera, Crosso- 
dactylus, Hylodes, Megaelosia, and Phantasmarana, with 
the latter being proposed here to accommodate six species 
assigned previously to the genus Megaelosia, which is now 
a monotypic genus. 

The inclusion in our phylogenetic analyses of six of the 
seven species previously assigned to Megaelosia allowed 
us to revisit the systematics of this genus and elucidate the 
interspecific relationships of M. goeldii. Phylogenetic and 
morphological evidence supports the definition of M. goel- 
dii as the only species of the genus Megaelosia, while the 
remaining species assigned historically to this genus in fact 
form a distinct clade. The monotypic Megaelosia is now 
demonstrated to be a sister group of the genus Hylodes, as 
inferred from our two phylogenetic analyses. 

In addition to molecular evidence, a number of mor- 
phological characters distinguish Phantasmarana from 
Hylodes and Megaelosia. As proposed by NUIN & VAL 
(2005), two synapomorphies of Hylodes may be recognized 
- the fusion of the prootic and the frontoparietals, and the 
rudimentary alary process of the hyoid. Megaelosia goeldii 
has two autapomorphic conditions, as described by Extas- 
Costa et al. (2017), that is, the absence of a vocal sac and 
the apical supplementary elements of the m. intermand- 
ibularis. Furthermore, NUIN & VAL (2005) described the 
zygomatic ramus as being larger than the optic ramus in 
M. goeldii, a trait that is absent in Hylodes. Phantasmarana 
bocainensis (missing from our molecular analysis) was in- 
cluded in the new genus due to the fact that the fold of the 
fifth toe does not reach the outer metatarsal tubercle (see 
GIARETTA et al. 1993), one of the putative synapomorphies 
of Phantasmarana. 

From a biogeographic perspective, the species of Phan- 
tasmarana occur in the Serra da Mantiqueira and Serra do 
Mar mountain ranges (in SP, RJ, and MG). The Serra da 
Mantiqueira hosts several endemic anuran taxa (CRUZ & 
FEIO 2007, Happap et al. 2013, SILVA et al. 2018), and phylo- 
geographic evidence supports the historical role of geomor- 
phological changes as a source of the vicariant processes 
that likely resulted in the high level of interspecific diver- 
sification observed in this region (CHavEs et al. 2014). The 
Phantasmarana populations, which are associated with high 
altitude, appear to represent a promising model for phyloge- 
ographic studies, which may further corroborate evidence 
of ongoing evolutionary processes in eastern Brazil. 
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No evidence of exclusive morphological characteris- 
tics of the tadpoles or adults have been identified in either 
of the well-supported clades recovered in Phantasmarana 
(P. massarti + P. boticariana) and (P. jordanensis + (P. apua- 
na + P. lutzae)) in previous studies (GIARETTA et al. 1993, 
POMBAL et al. 2003, NOGUEIRA-CosTa et al. 2012, SICHIERI 
et al. 2020). Although P massarti and P. boticariana occur 
at geographically proximal sites, the phylogeographic rela- 
tionship between them is unclear, given that they are found 
in distinct geological formations, with P massarti occur- 
ring in the Serra do Mar range (GIARETTA et al. 1993, FROST 
2021) and P. boticariana in the Serra da Mantiqueira (Mus- 
CAT et al. 2020a). The phylogeographic relationships of the 
clade composed of P. jordanensis, P. lutzae, and P. apuana 
are also unclear, although P. jordanensis and P. lutzae are 
both found in the Serra da Mantiqueira (NOGUEIRA-Cos- 
TA et al. 2012, IZECKSOHN & GOUVEA 1987), while P. apua- 
na occurs in the Serra do Caparaó (an extension of the Ser- 
ra da Mantiqueira, but approximately 460 km from the to- 
potype of P. jordanensis and 360 km from that of P. lutzae; 
SANTOS et al. 2011, ZORNOSA-ToRRES et al. 2020). Further 
research is needed to elucidate the phylogenetic relation- 
ships between P. bocainensis and the other Phantasmarana 
species. 

The diploid number of 2n = 26 is a putative plesiomor- 
phic condition of the family Hylodidae, as it is found in 
the genera Crossodactylus, Hylodes, and Megaelosia. In the 
most parsimonious scenario, the most recent common an- 
cestor of Phantasmarana would have had 2n = 28 chromo- 
somes, derived from a rearrangement of the 2n = 26 karyo- 
type. Subsequent rearrangements would have led to further 
increases in diploid numbers, as observed in P. boticariana, 
with 2n = 30 chromosomes (GIARETTA & AGUIAR-JR. 1998, 
Rosa et al. 2003) and, in the other clade, P. jordanensis, 
with 2n = 32 (Rosa et al. 2003). However, a more conclu- 
sive interpretation of the evolutionary changes in Phantas- 
marana would require chromosomal data from P. apuana, 
P lutzae, and P. bocainensis. 

The findings of the present study provide a better un- 
derstanding of the phylogenetic relationships in the Hylo- 
didae, in addition to unravelling the historical taxonomic 
uncertainties of Megaelosia, with most species now being 
assigned to the new genus, Phantasmarana. The present 
study thus contributes to a more comprehensive viewpoint 
on the evolutionary processes that have formed this fas- 
cinating and diverse anuran family, and it will provide an 
essential database for future studies in the mountains of 
eastern Brazil. 


Acknowledgements 


We thank L. Tactou1, C. LAMBERTINI, C. ZORNOSA-TORRES, 
B. Dias, L. RIBEIRO, J. RUGGERI, J. PRADO, J. BOVOLON, and 
A. BERTOLAZZzI for their assistance during fieldwork, and K. 
SEGER, R. TENORIO, A. F. pa Costa, R. N. PEDROSO, and K. P. 
GATTO for their extensive help in the laboratory. We are thankful 
for J. KOHLER and B. RöppeR for leading the excellent review- 
ing and editorial process. We are also grateful to the administra- 


282 


tion of Caparaó National Park for granting access to the park, and 
their logistic support of our research. GAA and LFT are grateful 
to the São Paulo Research Foundation (FAPESP #2016/25358-3, 
#2019/03170-0) and the Brazilian National Council for Scientific 
and Technological Development (CNPq #300896/2016-6). 


References 


AGUIAR-JR., O., K. A. CARVALHO, A. A. GIARETTA & S. M. REC- 
CO-PIMENTEL (2003): Cytogenetics of Hylodes and Crossodac- 
tylus species (Anura, Leptodactylidae) with comments on Hy- 
lodinae/Dendrobatid relationships. - Genetica, 121: 43-53. 


ALVES, G. A & L. F ToLepo (2017): Food items of the rheophilic 
Torrent Frog Megaelosia apuana Pombal et al., 2003 (Anura, 
Hylodidae). - Herpetology Notes, 10: 535-536. 


AMARO-GHILARDI, R. C. (2005): Estudos citogenéticos compara- 
tivos em 26 espécies da familia Leptodactylidae (Amphibia, 
Anura). - Unpubl. Ph.D. thesis. 


AucGusTo-ALVES, G., S. A. DENA & L. F. TOLEDO (2018): Visual 
communication and aggressive behaviour in a giant mute tor- 
rent-frog, Megaelosia apuana (Anura; Hylodidae). - Amphi- 
bia-Reptilia, 39: 260-264. 

BAUMANN, F. (1912): Brasilianische Batrachier des Berner 
Naturhistorischen Museums nebst Untersuchungen über die 
geographische Verbreitung der Batrachier in Brasilien. Zoo- 
logische Jahrbücher. - Abteilung fiir Systematik, Geographie 
und Biologie der Tiere, 33: 87-172. 


BEÇAK, M. L. (1968): Chromosomal analysis of eighteen species 
of Anura. - Caryologia 21: 191-208. 


BLOTTO, B. L., J. J. NuÑez, N. G. Basso, C. A. ÚBEDA, W. C. 
WHEELER & J. FAIVOVICH (2012): Phylogenetic relationships 
of a Patagonian frog radiation, the Alsodes + Eupsophus clade 
(Anura: Alsodidae), with comments on the supposed para- 
phyly of Eupsophus. - Cladistics, 29: 113-131. 

BOGART, J. P. (1970): Systematic problems in the amphibian Fam- 
ily Leptodactylidae (Anura) as indicated by karyotypic analy- 
sis. - Cytogenetics, 9: 369-383. 

BOGART, J. P. (1991): The influence of life history on karyotypic 
evolution in frogs. - pp. 233-258 in: GREEN D. M. & S. K. SEs- 
sons (eds): Amphibians cytogenetics and evolution. - Aca- 
demic Press, San Diego. 


BOUCKAERT, R., T. G. VAUGHAN, J. BARIDO-SOTTANI, S. DU- 
CHÊNE, M. FOURMENT, A. GAVRYUSHKINA, J. HELED, G. JONES, 
D. KUHNERT, N. DE Mato, M. MATSCHINER, F. K. MENDES, 
N. F. MÜLLER, H. A. OGILVIE, L. pu PLEssIs, A. POPINGA, 
A. RAMBAUT, D. RASMUSSEN, I. SIVERONI, M. A. SUCHARD, 
C-H. Wu, D. XIE, C. ZHANG, T. STADLER, & A. J. DRUMMOND 
(2019): BEAST 2.5: An advanced software platform for Baye- 
sian evolutionary analysis. - PLOS Computational Biology, 
15(4): €1006650. 


BRUM-ZORRILLA, N & F. A. SAEZ (1968): Chromosomes of Lepto- 
dactylidae (Amphibia—Anura). - Experientia, 27: 969. 

CARVALHO, T. R., K. R. SEGER, F. M. MAGALHÃES, L. B. LOURENCO 
& C. F. B. HaAppap (2020): Systematics and cryptic diversifica- 
tion of Leptodactylus frogs in the Brazilian campo rupestre. - 
Zoologica Scripta, oo: 1-18. 


CHAVES, A. V., G. H. S. Freitas, M. F. VASCONCELOS, & F. R. SAN- 
TOS (2014): Biogeographic patterns, origin and speciation of 
the endemic birds from eastern Brazilian mountaintops: a re- 
view. - Systematics and Biodiversity, 13: 1-16. 


Paraphyly in the giant torrent-frogs and description of a new genus 


Cruz, C. A. G. & R. N. FEIo (2007): Endemismos em anfíbios em 
areas de altitude na Mata Atlantica no sudeste do Brasil. - pp. 
117-126, in: NASCIMENTO, L. B. & M. E. OLIVEIRA (eds): Her- 
petologia no Brasil II. - Sociedade Brasileira de Herpetologia, 
Belo Horizonte. 


DE Lucca, E. J., J. JIM & F. Foresti (1974): Chromosomal studies 
in twelve species of Leptodactylidae and one Brachycephali- 
dae. - Caryologia, 27: 183-191. 

DE WITTE, G. F (1930): Liste des reptiles et batraciens récoltés 
au Brésil par la Mission Massart (1922-23) et description de 
sept nouvelles espèces. - pp. 213-230 in: Massarr, J. (ed.): Une 
Mission Biologique Belge au Brésil (août 1922—mai 1923) par 
Jean Massart, Raymond Bouillene, Paul Ledoux, Paul Brien, 
et Albert Navez. Volume 2. - Impremerie Médicale et Scienti- 
fique, Bruxelles. 


DENARO, L. (1972): Karyotypes of Leptodactylidae anurans. - 
Journal of Herpetology, 6: 71-74. 


Dusots, A., A. OHLER & R. A. PYRON (2021): New concepts and 
methods for phylogenetic taxonomy and nomenclature in zo- 
ology, exemplified by a new ranked cladonomy of recent am- 
phibians (Lissamphibia). - Megataxa, 5: 1-738. 

DumÉRIL, A. M. C. & L. BIBRON (1841): Erpétologie Genérale ou 
Histoire Naturelle Complète des Reptiles. Volume 8. - Librarie 
Enclyclopedique de Roret, Paris. 


Exras-Costa, A. J., R. MONTESINOS, T. GRANT & J. FAIVOVICH 
(2017): The vocal sac of Hylodidae (Amphibia, Anura): Phylo- 
genetic and functional implications of a unique morphology. 
- Journal of Morphology, 278: 1506-1516. 


Faivovicn, J., C. F. B. Happap, P. C. A. Garcia, D. R. FROST, 
J. A. CAMPBELL & W. C. WHEELER (2005): Systematic review 
of the frog family Hylidae, with special reference to Hylinae: 
phylogenetic analysis and taxonomic revision. - Bulletin of 
the American Museum of Natural History, 294: 1-240. 


FELLER, A & S. B. HEDGEs (1998): Molecular evidence for the ear- 
ly history of living amphibians. - Molecular Phylogenetics and 
Evolution, 9: 509-516. 


FELSENSTEIN, J. (2004): Inferring phylogenies. - Sinauer Associ- 
ates, Sunderland. 


FITZINGER, L. J. F. J. (1826): Neue Classification der Reptilien 
nach ihren Natiirlichen Verwandtschaften nebst einer Ver- 
wandtschafts-Tafel und einem Verzeichnisse der Reptilien- 
Sammlung des K. K. Zoologischen Museums zu Wien. - J. G. 
Heubner, Wien. 


Frost, D. R. (2021): Amphibian species of the world: an online 
reference. Version 6.1. Electronic database accessible at Amer- 
ican Museum of Natural History. - American Museum of 
Natural History, New York. Available at: https://amphibian- 
softheworld.amnh.org/index.php, accessed 8 January 2021. 


GATTO, K. P., J. V. MATTOS, K. R. SEGER & L. B. LOURENÇO (2018): 
Sex Chromosome Differentiation in the Frog Genus Pseudis 
Involves Satellite DNA and Chromosome Rearrangements. - 
Frontiers in Genetics, 9: 301. 


GIARETTA, A. A., W. C. A. BOKERMAN & C. F. B. HADDAD (1993): 
A review of the genus Megaelosia (Anura: Leptodactylidae) 
with a description of a new species. - Journal of Herpetolo- 
gy» 27: 276-285. 

GLICK, L. & I. MayROSE (2014): ChromEvol: assessing the pattern 
of chromosome number evolution and the inference of poly- 
ploidy along a phylogeny. - Molecular Biology and Evolution, 
31.7: 1914-1922. 


GIARETTA, A. A. & O. AGUIAR-JR. (1998): A new species of Mega- 
elosia from the Mantiqueira Range, Southeastern Brazil. - 
Journal of Herpetology, 32: 80-83. 


GOLOBOFF, P. A. & S. A. CATALANO (2016): TNT version 1.5, in- 
cluding a full implementation of phylogenetic morphomet- 
rics. — Cladistics, 32: 221-238. 


GRANT, T., M. Rapa, M. ANGANOY-CRIOLLO, A. BATISTA, P. H. 
Dias, A. M. JECKEL, D. J. MACHADO & J. V. RUEDA-ALMONA- 
CID (2017): Phylogenetic systematics of dart-poison frogs and 
their relatives revisited (Anura: Dendrobatoidea). - South 
American Journal of Herpetology, 12(s1): S1-Sgo. 


GRAYBEAL, A. (1997): Phylogenetic relationships of bufonid frogs 
and tests of alternate macroevolutionary hypotheses charac- 
terizing their radiation. - Zoological Journal of the Linnean 
Society, 119: 297-338. 


GREEN, D. M & S. K. Sessions (1991): Nomenclature for Chromo- 
somes. — pp. 431-432 in: GREEN, D. M & S. K. Sessions (eds): 
Amphibian cytogenetics and evolution. - Academic Press, San 
Diego. 

HADDAD, C. F B., L. E TOLEDO, C. P. A. PRADO, D. LOEBMANN, 
J. L. GASPARINI & I. SAZIMA (2013): Guia dos anfíbios da Mata 
Atlantica - Diversidade e Biologia. - Anolis Books Editora, 
São Paulo. 


HALL, T. A. (1999): BioEdit: a user-friendly biological sequence 
alignment editor and analysis program for Windows 95/98/ 
NT. - Nucleic Acids Symposium Series, 41: 95-98. 


HEDGES, S. B. (1994): Molecular evidence for the origin of birds. 
- Proceedings of the National Academy of Sciences, USA, 91: 
2621-2624. 


HEYER, W. R. (1983): Variation and systematics of frogs of the ge- 
nus Cycloramphus (Amphibia, Leptodactylidae). - Arquivos 
de Zoologia, 30: 235-339. 


HowELL, W. M. & D. A. BLACK (1980): Controlled silver staining 
of nucleolus organizer regions with a protective colloidal de- 
veloper: a 1-step method. - Experientia, 36: 1014-1015. 


HUELSENBECK, J. P. & D. M. HILtis (1993): Success of phyloge- 
netic methods in the four-taxon case. — Systematic Biology, 
42: 247-264. 

IZECKSOHN, E. & E. GouvEA (1987 “1985”): Nova especie de Mega- 
elosia de Itatiaia, Estado do Rio de Janeiro. - Arquivos de Uni- 
versidade Federal Rural do Rio de Janeiro, 8: 17-22. 


Katou, K. J., J. ROZEWICKI & K. D. YAMADA (2019): MAFFT 
online service: multiple sequence alignment, interactive se- 
quence choice and visualization. - Briefings in Bioinformat- 
ics, 20(4): 1160-1166. 


Kina, M. (1980): C-banding studies on Australian hylid frogs: 
secondary constriction structure and the concept of euchro- 
matin transformation. - Chromosoma, 80: 191-217. 


KING, M & R. Rore (1976): Karyotypic variation in the australian 
gekko Phyllodactylus Marmoratus (Gray) (Gekkonidae, Rep- 
tilia). - Chromosoma, 54: 75-87. 

KUMAR, S., G. STECHER, M. Li, C. KNyAz & K. TAMURA (2018): 
MEGA X: Molecular evolutionary genetics analysis across 
computing platforms. - Molecular Biology and Evolution, 
35(6): 1547-1549. 

MELO, A. S., S. M. RECCO-PIMENTEL & A. A. GIARETTA (1995): 
The karyotype of the stream dwelling frog Megaelosia mas- 
sarti (Anura, Leptodactylidae, Hylodinae). - Cytologia, 60: 
49-52. 


283 


STENIO EDER VITTORAZZI et al. 


MIRANDA-RIBEIRO, A. (1923): Elosia Tsch. e os generos correlatos. 
— Revista do Museu Paulista. São Paulo, 13: 813-821. 


Muscat, E., D. STUGINSKI, I. Nunes, I. MARTINS, G. AuGus- 
TO-ALVES, S. E. VITTORAZZI, L. F ToLepo & M. T. Morott. 
(2020a): Update on the geographic distribution of three poor- 
ly known frog species in the Mantiqueira mountain range. - 
Herpetology Notes, 13: 573-577. 


Muscat, E., D. STUGINSKI, M. T. MororI & L. F. TOLEDO (2020b): 
An unidentified call emission by Megaelosia boticariana. - 
Herpetology Notes, 13: 805-808. 


NOGUEIRA-COSTA, P., P. ALMEIDA-SANTOS, C. A. G. CRUZ & U. 
CARAMASCHI (2012): The giant tadpoles of Megaelosia jor- 
danensis (Heyer, 1983). - Zootaxa, 3581(1): 86-88. 


Nurin, P. A. S. & F. C. VAL (2005): Phylogenetic analysis of the 
subfamily Hylodinae (Anura, Leptodactylidae) based on mor- 
phological characters. - Amphibia-Reptilia, 26: 139-147. 


NYLANDER, J. A. A. (2004): MrModeltest 2.2. Program distrib- 
uted by the author. - Evolutionary Biology Centre, Uppsala 
University, Upsalla. 


PeLoso, P. L. V. & S. E. PAvAN (2007): Natural history notes. 
Megaelosia apuana. Prey. - Herpetological Review, 38: 188. 


POMBAL, J. P., G. M. PRADO & C. CANEDO (2003): A new species 
of giant torrent frog, genus Megaelosia, from the Atlantic Rain 
Forest of Espirito Santo, Brazil (Amphibia: Leptodactylidae). - 
Journal of Herpetology, 37: 453-460. 


PaLumBI, S. R., A. Martin, W. O. McMILLan, L. STICE & G. 
GRABOWSKI (2002): The simple fool’s guide to PCR, Version 
2.0. — University of Hawaii, 45 pp. 


Pyron, R. A & J. J. WIENS (2011): A large-scale phylogeny of Am- 
phibia including over 2,800 species, and a revised classifica- 
tion of extant frogs, salamanders, and caecilians. - Molecular 
Phylogenetics and Evolution. 61: 543-583. 


QUERCIA, C. A., E. Y. SUAREZ-VILLOTA, F. FORESTI & J. J. NUNEZ 
(2020): Comparative cytogenetics of the ground frogs Eu- 
psophus emiliopugini Formas, 1989 and E. vertebralis Grandi- 
son, 1961 (Alsodidae) with comments on their inter- and in- 
traspecific chromosome differentiation. - Comparative Cy- 
togenetics, 14(1): 61-74. 


RAMBAUT, A., A. J. DRUMMOND, D. XIE, G. BAELE & M. A. Su- 
CHARD (2018): Posterior summarisation in Bayesian phylo- 
genetics using Tracer 1.7. - Systematic Biology, 67: 901-904. 
Available at: http://beast.community/tracer. 


RAMBAUT, A. (2012): FigTree v1.4.2. - Available at: http://tree.bio. 
ed.ac.uk/software/figtree/, accessed 11 September 2020. 


RICE, A. & I. MAyRosE (2020): Model Adequacy Tests for Prob- 
abilistic Models of Chromosome-Number Evolution. - New 
Phytologist. 

Rosa, C., O. AGUIAR-JR., A. A. GIARETTA & S. M. RECCO-PIMEN- 
TEL (2003): Karyotypic variation in the genus Megaelosia (An- 
ura, Leptodactylidae) with the first description of a B chromo- 
some in leptodactylid frogs. - Copeia, 2003(1): 166-174. 

SANTOS, P. A., E. T. Sttva, B. H. B. FELHBERG, M. T. T. SANTOS & 
P. C. A. Garcia (2011): Amphibia, Anura, Hylodidae, Mega- 
elosia apuana Pombal, Prado and Canedo, 2003: Distribution 
extension, new state record and geographic distribution map. 
- Check List, 7(4): 394-396. 


SEGALLA, M. V., B. BERNECK, C. CANEDO, U. CARAMASCHI, A. G. 
Cruz, P. C. A. GARCIA, T. GRANT, C. F. B. HADDAD, A. C. C. 
LOURENÇO, S. MANGIA, T. Mott, F. P. WERNECK, L. B. NASCI- 


284 


MENTO, L. F. TOLEDO & J. A. LAGONE (2021): List of Brazilian 
amphibians. - Herpetologia Brasileira, 10(1): 118-208. 


SICHIERI, G., M. T, Moroti, F. R. Costa, E. MUSCAT & I. NUNES 
(2020): Redescription of the tadpole of rare giant torrent frog 
Megaelosia boticariana Giaretta & Aguiar, 1998 (Anura, Hylo- 
didae) with comments on natural history. - Zootaxa, 4896 (1): 
140-144. 

SILVA, E. T., M. A. A. PEIXOTO, F. S. F. LEITE, R. N. FEIO & P. C. 
A. GARCIA (2018): Anuran distribution in a highly diverse re- 
gion of the Atlantic Forest: The Mantiqueira Mountain Range 
in southeastern Brazil. - Herpetologica, 74 (4): 294-305. 


Titus, T. A. (1992): A Phylogenetic Analysis of the Desmognathi- 
nae (Caudata: 1229 Plethodontidae): Evolutionary Patterns In- 
ferred from Mitochondrial DNA 1230 Sequences. - Ph.D. the- 
sis., Univ. of Kansas at Lawrence. 


TOLEDO, L. F, C. G. BECKER, C. F. B. HADDAD & K. R. ZAMUDIO 
(2014): Rarity as an indicator of endangerment in Neotropical 
frogs. - Biological Conservation, 179: 54-62. 


WuIrTeE, M. J. D. (1954): Animal Cytology and Evolution. - Cam- 
bridge University Press, Cambridge. 


ZORNOSA-TORRES C., G. AUGUSTO-ALVES, M. L. Lyra, J. C. DA 
S. JÚNIOR, P. C. A. GARCIA, F. LEITE, V. VERDADE, M. T. Ro- 
DRIGUES, J. L. GASPARINI, C. F. B. HADDAD & L. F TOLEDO 
(2020): Anurans of the Caparaó National Park and surround- 
ings, southeast Brazil. - Biota Neotropica, 20(3): e20190882. 


Supplementary data 


The following data are available online: 


Supplementary Table S1. Specimens included in our analysis with 
their localities and identification (ID) numbers in scientific col- 
lections and the GenBank accession numbers of the nucleotide 
sequences, with respective bibliographic citations. 


Supplementary Table $2. Measurements of analyzed individuals, 
including information available on literature. Abbreviations: SVL 
= snout-vent length; TL = total length; HW = head width; HL = 
head length; FDW = finger discs width. 

Appendix S1. Additional hylodid specimens or recordings exam- 


ined, to complement the information available in the literature, 
which are deposited in the following Brazilian collections. 


